home *** CD-ROM | disk | FTP | other *** search
Wrap
Frequently Asked Questions (FAQS);faqs.250 General index for the MATLAB User Group software library ---------------------------------------------------------- Currently there are the following subdirectories: approximation approximation theory archive old MATLAB user group digests control control theory dataanalysis data analysis and statistics graphics graphics programs integration numerical integration linearalgebra linear algebra utilities misc miscellaneous ode ordinary differential equations optimization as the name says pde partial differential equations rootfinding zero-finding routines specialfunctions special functions teaching for classroom use tools miscellaneous tools In order to get an index for a subdirectory (tools, say) send the message send index from MATLAB/tools to netlib@ornl.gov. In order to get some code, (unbundle in the `tools' directory, say), send the message send unbundle from MATLAB/tools to netlib@ornl.gov. FOR STUDENTS: Prentice Hall has published a student edition of matlab which contains a book and set of disks for PCs and Macs. The software is limited only in matrix size (32 x 32 matrix; 1024 elements) and in its ability to import or call C or Fortran subroutines. On the plus side, it is able to run without a coprocessor (it will use one if it is present) and it includes a subset of the Signal Processing and Controls Toolboxes, The Signals and Systems Toolbox, which provides for added functionality. Book only (about US$30): ISBN = 0138560064; Book + disk: (about US$50 ISBN=0-13-855974-0 for 3.5" or ISBN=0-13-855982-1 for 5.25 Macintosh version: ISBN=0-13-855990-2. There will be related books out by mid to late 1993 : Computer Aided Signal Processing with MATLAB, by Burrus, Oppenheim, McClellan, Parks, Schafer, and Schussler; and Signal Processing : A Computer Approach, by Etter. More books in this MATLAB Curriculum Series are planned. For general info: matlab@prenhall.com [From the Matlab Users Group (Editor, hwilson@ua1vm.ua.edu)] -------------------------------------------------- Q1.2.6: What are the Signal Processing Packages for Mathematica? Where can I get them? Package: Signal Processing Pacakages (SPP) and Notebooks. Description: Public domain extensions to Mathematica. Enables the symbolic manipulation of signal processing expressions: 1-D discrete/continuous convolutions and 1-D/m-D linear transforms (Laplace, Fourier, z, DTFT, and DFT). For linear transforms, you can specify your own transform pairs and see the intermediate computations. Great for showing students how to take transforms, or for deriving input-output relationships in a transform domain. Additional abilities include analog filter design, solving DE's using transforms, converting signal processing expressions to their equivalent TeX forms, and number theoretic operations (Bezout numbers, Smith Form decompositions, and matrix factors). Accompanying SPP are tutorial notebooks on analog filter design, Fourier analysis, piecewise convolution, and the z-transform (includes a discussion of fundamentals of digital filter design). These Notebooks illustrate difficult concepts (such as the flip-and-slide view of convolution) through animation. A Notebook reader is available in the public domain for Macintosh computers (a BinHexed version is on the ftp site). FOR STUDENTS: A student version of Mathematica is available for $175. The price includes a copy of the reference manual. The only drawbacks to the student version are that the floating point coprocessor is dis- abled and that upgrades cannot be ordered. Contact: Brian Evans, evans@eedsp.gatech.edu. To obtain: anonymous ftp to gauss.eedsp.gatech.edu (130.207.226.24). (available in UNIX, Mac, and MS-DOS archive formats). -------------------------------------------------- Q1.2.7: What is the Control Systems Analysis Packages for Mathematica? Where can I get them? Package: Control Systems Analysis Package (COSYPAK) and Notebooks Description: Public domain extension to Mathematica. Classical and state-space control analysis and design methods. The Notebooks supplement the material in the textbook "Modern Controls Theory" by Ogata. Largely based on the Signal Processing Packages (SPP, see above). Contact: Dr. Sreenath, sree@veda.esys.cwru.edu. To obtain: anonymous ftp veda.esys.cwru.edu (129.22.40.9). -------------------------------------------------- Q1.2.8: What are some other Mathematica DSP Notebooks? The following Mathematica notebooks can be ftped from ccrma-ftp.stanford.edu: pub/DSP/GenHamming.ma.Z Generalized Hamming windows pub/DSP/Kaiser.ma.Z The Kaiser window pub/DSP/WinFlt.ma.Z Digital filter design by the "window method" (There are other DSP related items in pub/DSP on ccrma-ftp; see other sections of this FAQ for details). -------------------------------------------------- Q1.2.9: What is the Linear Systems Toolbox for Maple? Where can I get it? Package: Linear systems toolbox for Maple. Description: Public domain extension to Maple. Contact: Tony Richardson, amr@mpl.ucsd.edu. To obtain: anonymous ftp to cameron.egr.duke.edu, file pub/maple/linsys1.2.tar.Z. -------------------------------------------------- Q1.2.10: Where can I get text to speech conversion software? Free (but not public domain) text to speech conversion software is available via anonymous ftp from wilma.cs.brown.edu in the pub directory as speak.tar.Z. It will compile and run on a SPARC's built-in audio after modifying speak.c with the path of your libaudio.h (e.g., /usr/demo/SOUND/libaudio.h). It's a simple phoneme concatenation system with commensurate synthesized speech quality (a directory of phoneme audio files is included). [Joe Campbell, jpcampb@afterlife.ncsc.mil] -------------------------------------------------- Q1.2.11: Where can I get filter design software? One source is the following: in the August 92 issue of IEEE Transactions on Signal Processing there is a paper entitled "METEOR: A Constraint-Based FIR Filter Design Program" by Kenneth Steiglitz, Thomas W. Parks and James F. Kaiser. They describe an FIR design program which allows specification of the target frequency response characteristics in a fairly generalised and flexible way. As well as designing filters, the program can optimise filter lengths and push band limits. The paper contains a footnote which says "Pascal and C versions of source code are available to anonymous ftp at princeton.edu in the directory /pub as meteor.p, form.p, meteor.c and form.c". True, they are. They appear to work. The Pascal versions have been put through p2c to get the C versions; all the needed Pascal library stuff is included in the C versions and they built error-free out of the box for me on an SGI machine. One catch is, there is no manual - you need the paper to know how to drive the programs. [Steve Clift, clift@ml.csiro.au] { There are other free filter design programs floating around out there, such as optfir/wfir. Does anyone know of ftp sites? } ============================================================================= 2. Algorithms and standards. This section deals with DSP algorithms and related standards. ======================================================================= Q2.1: Where can I get some algorithms for general DSP? The following archives contain things such as matrix operations, FFT's and generally useful things like that, as opposed to complete applications: (1) host ux1.cso.uiuc.edu - log in as anonymous, get the file math/README. It summarises the contents of the archive. (2) Netlib, which serves some of this software via email. Try mail to netlib@ORNL.GOV with "send help" in the subject field. For Europe: Internet: netlib@nac.no EARN/BITNET: netlib%nac.no@norunix.bitnet X.400: s=netlib; o=nac; c=no; EUNET/uucp: nac!netlib For the Pacific, try netlib@draci.cs.uow.edu.au For background about netlib, see Jack J. Dongarra and Eric Grosse, "Distribution of Mathematical Software Via Electronic Mail," Comm. ACM (1987) 30,403--407. A similar collection of statistical software is available from statlib@temper.stat.cmu.edu. The symbolic algebra system REDUCE is supported by reduce-netlib@rand.org. (3) The Naval Surface Warfare Center has a library of mathematical Fortran subroutines that may be of use. From the report itself: NSWC Library of Mathematical Subroutines Report No.: NSWC TR 90-21, January 1990 by Alfred H. Morris, Jr. Naval Surface Warfare Center (E43) Dahlgren, VA 22448-5000 U.S.A. Distribution: Approved for public release; distribution unlimited. Abstract: The NSWC library is a library of general-purpose Fortran subroutines that provide a basic computational capability in a variety of mathematical activities. Emphasis has been placed on the transportability of the codes. Subroutines are available in the following areas: Elementary Operations, Geometry, Special Functions, Polynomials, Vectors, Matrices, Large Dense Systems of Linear Equations, Banded Matrices, Sparse Matrices, Eigenvalues and Eigenvectors, l1 Solution of Linear Equations, Least-Squares Solution of Linear Equations, Optimization, Transforms, Approximation of Functions, Curve Fitting, Surface Fitting, Manifold Fitting, Numerical Integration, Integral Equations, Ordinary Differential Equations, Partial Differential Equations [Witold Waldman, witold@hotblk.aed.dsto.gov.au] { Does anyone know an FTP site for this library? } If you don't know where to find what you're after, try archie. ======================================================================= Q2.2: What are CELP and LPC? Where can I get the source for CELP and LPC? CELP stands for "code excited linear prediction". LPC stands for "linear predictive coding". They are compression algorithms used for low bit rate (2400 and 4800 bps) speech coding. The U.S. DoD's Federal Standard 1016 (FS 1016) based 4800 bps code excited linear prediction voice coder version 3.2 (CELP 3.2) Fortran and C simulation source codes are available for worldwide distribution at no charge (on DOS diskettes, but configured to compile on Sun SPARC stations) from: Bob Fenichel National Communications System Washington, D.C. 20305 1-703-692-2124 1-703-746-4960 (fax) Example input and processed speech files, a technical information bulletin, and the official standard "Federal Standard 1016, Telecommunications: Analog to Digital Conversion of Radio Voice by 4,800 bit/second Code Excited Linear Prediction (CELP)" are included at no charge. Unfortunately, a document that is a vital part of the CELP release package is not available in electronic form. Anyone serious interested in CELP should obtain the document: Details to Assist in Implementation of Federal Standard 1016 CELP. National Communications System, Office of Technology & Standards, 1992. Technical Information Bulletin 92-1. It is available from the above address. The 4800 bps CELP code is available via anonymous FTP on furmint.nectar.cs.cmu.edu (128.2.209.1111) in celp.audio.compression. It is unsupported, and has not been compiled or tested. This version takes a lot of CPU time. It is about 15 times slower than real-time. That is, to code 1 minute of speech takes about 15 minutes. It has directions on how to reduce the codebook search to make it faster, but they are only for the speech coding expert. The following article describes the FS 1016 4.8-kbps CELP coder: Campbell, Joseph P. Jr., Thomas E. Tremain and Vanoy C. Welch, "The Proposed Federal Standard 1016 4800 bps Voice Coder: CELP," Speech Technology Magazine, April/May 1990, p. 58-64. The U.S. DoD's Federal-Standard-1015/NATO-STANAG-4198 based 2400 bps linear prediction coder version 53 (LPC-10e v53) Fortran or C simulation source codes are available on a limited basis upon written request to: Tom Tremain Department of Defense Ft. Meade, MD 20755-6000 USA The U.S. Federal Standard 1015 (NATO STANAG 4198) is described in: Thomas E. Tremain, "The Government Standard Linear Predictive Coding Algorithm: LPC-10," Speech Technology Magazine, April 1982, p. 40-49. Copies of the official standards FS 1015 and FS 1016 are available for US $2.50 each from: GSA Rm 6654 7th & D St SW Washington, D.C. 20407 1 (202) 708-9205 Realtime DSP code for FS-1015 and FS-1016 is sold by several vendors, including: John DellaMorte DSP Software Engineering 165 Middlesex Tpk, Suite 206 Bedford, MA 01730 1-617-275-3733 1-617-275-4323 (fax) dspse.bedford@channel1.com and Analogic Systems 2916 Ramona Palo Alto, CA 94306 (415) 323-3232 DSP Software Engineering's FS-1016 code can run on a DSP Research's Tiger 30 or on Intellibit's AE2000 TMS320C31 based 3" by 2.5" card. See section 4.1 for more on these cards. Analogic's product runs on a 27 MHz DSP56001 chip. [Most of the above from Joe Campbell, jpcampb@afterlife.ncsc.mil, with additions from DanFrankowski, drankow@cs.umn.edu, and Ed Hall, edhall@rand.org] ======================================================================= Q2.5: What is ADPCM? Where can I get source for it? ADPCM stands for Adaptive Delta Pulse Code Modulation. It is a standard for speech compression and decompression. It takes 16-bit linear PCM samples samples and converts them to 4-bit samples, so the compression rate is 4x. There is public domain C code available via anonymous ftp at ftp.cwi.nl in /pub/adpcm.shar written by Jack Jansen (email Jack.Jansen@cwi.nl). It is very programmer-friendly. The ADPCM code used is the Intel/DVI ADPCM code which is being recommended by the IMA Digital Audio Technical Working Group. It allows the following calls: adpcm_coder(short inbuf[], char outbuf[], int nsample, struct adpcm_state *state); adpcm_decoder(char inbuf[], short outbuf[], int nsample, struct adpcm_state *state); The routines have been tested on an SGI Indigo running Irix 4.0.2 and on a Sparcstation 1+ running SunOS 4.1.1. On a Sun, the code will compress at 250Ksample/sec and decompress at 300Ksample/sec. On an SGI, the compressor runs at 350Ksample/sec and the decompressor at 700Ksample/sec. Note that this is NOT a CCITT G722 coder. The CCITT ADPCM standard is much more complicated, probably resulting in better quality sound but also in much more computational overhead. [From Dan Frankowski, drankow@cs.umn.edu; Jack Jansen, Jack.Jansen@cwi.nl] ======================================================================= Q2.4: How does pitch perception work, and how do I implement it on my DSP chip? Pitch is officially defined as "That attribute of auditory sensation in terms of which sounds may be ordered on a musical scale." Several good examples illustrating the subtleties of pitch perception are included in the "Auditory Demonstrations CD" which is available from the Acoustical Society of America, Woodbury, NY 10797 for $20. A good general reference about the psychology of pitch perception is the book: B.C.J. Moore, "An Introduction to the Psychology of Hearing", Academic Press, London, 1989. This book is available in paperback and makes a good desk reference. An algorithm implementation that matches a large body of psychoacoustical work, but which is computationally very intensive, is presented in the paper: Malcolm Slaney and Richard Lyon, "A Perceptual Pitch Detector," Proceedings of the International Conference of Acoustics, Speech, and Signal Processing, 1990, Albuquerque, New Mexico. The definitive papers describing the use of such a perceptual pitch detector as applied to the classical pitch literature is in: Ray Meddis and M. J. Hewitt. "Virtual pitch and phase sensitivity of a computer model of the auditory periphery. " Journal of the Acoustical Society of America 89 (6 1991): 2866-2682. and 2883-2894. The current work that argues for a pure spectral method starts with the work of Goldstein: J. Goldstein, "An optimum processor theory for the central formation of the pitch of complex tones," Journal of the Acoustical Society of America 54, 1496-1516, 1973. Two approaches are worth considering if something approximating pitch is appropriate. The people at IRCAM have proposed a harmonic analysis approach that can be implemented on a DSP Boris Doval and Xavier Rodet, "Estimation of Fundamental Frequency of Musical Sound Signals," Proceedings of the 1991 International Conference on Acoustics, Speech, and Signal Processing, Toronto, Volume 5, pp. 3657-3660. The classic paper for time domain (peak picking) pitch algorithms is: B. Gold and L. Rabiner, "Parallel processing techniques for estimating pitch periods of speech in the time domain," Journal of the Acoustical Society of America, 46, pp 441-448, 1969. Finally, a word of caution: Pitch is not single-valued. We can hear a sound and match it to several different pitches. Imagine the number of instruments in an orchestra, each with its own pitch. Even a single sound can have more than one pitch. See for example Demonstration 27 from the ASA Auditory Demonstrations CD. [The above from Malcolm Slaney, Apple Computer, and John Lazzaro, U.C. Berkeley.] ======================================================================= Q2.5: What standards are there for digital audio? What is AES/EBU? What is S/P-DIF? The "AES/EBU" (Audio Engineering Society / European Broadcast Union) digital audio standard is probably the most popular digital audio standard today. Most consumer and professional digital audio devices (CD players, DAT decks, etc.) that feature digital audio I/O support AES/EBU. AES/EBU is a bit-serial communications protocol for transmitting digital audio data through a single transmission line. It provides two channels of audio data (up to 24 bits per sample), a method for communication control and status information ("channel status bits"), and some error detection capabilities. Clocking information (i.e., sample rate) is derived from the AES/EBU bit stream, and is thus controlled by the transmitter. The standard mandates use of 32 kHz, 44.1 kHz, or 48 kHz sample rates, but some interfaces can be made to work at other sample rates. AES/EBU provides both "professional" and "consumer" modes. The big difference is in the format of the channel status bits mentioned above. The professional mode bits include alphanumeric channel origin and destination data, time of day codes, sample number codes, word length, and other goodies. The consumer mode bits have much less information, but do include information on copy protection (naturally). Additionally, the standard provides for "user data", which is a bit stream containing user-defined (i.e., manufacturer-defined) data. According to Tim Channon, "CD user data is almost raq CD subcode; DAT is StartID and SkipID. In progfessional mode, there is an SDLC protocol or, if DAT, it may be the same as consumer mode." There physical connection media are commonly used with AES/EBU: balanced (differential), using two wires and shield in three-wire microphone cable with XLR connectors; unbalanced (single-ended), using audio coax cable with RCA jacks; and optical (via fiber optics). "S/P-DIF" (Sony/Philips Digital Interface Format) typically refers to AES/EBU operated in consumer mode over unbalanced RCA cable. Note that S/P-DIF and AES/EBU mean different things depending on how much of a purist you are in the digital audio world; see the Finger article below. References: Finger, Robert, "AES3-199X: The Revised Two Channel Digital Audio Interface (DRAFT)", presented at the 91st Convention of the Audio Engineering Society, October 4-8, 1991. Reprints: AES, 60 East 42nd St., New York, NY, 10165. [The above from Phil Lapsley, phil@ohm.Berkeley.EDU, and Tim Channon, tchannon@black.demon.co.uk] ============================================================================= Q2.6: What is mu-law encoding? Where can I get source for it? Mu-law (also "u-law") encoding is a form of logarithmic quantization or companding. It's based on the observation that many signals are statistically more likely to be near a low signal level than a high signal level. Therefore, it makes more sense to have more quantization points near a low level than a high level. In a typical mu-law system, linear samples of 14 to 16 bits are companded to 8 bits. Most telephone quality codecs (including the Sparcstation's audio codec) use mu-law encoded samples. Desktop Sparc machines come with routines to convert between linear and mu-law samples. On a desktop Sparc, see the man page for audio_ulaw2linear in /usr/demo/SOUND/man. Craig Reese posted the source of similar routines to comp.dsp in August '92. References: CCITT Recommendation G.711 (very difficult to follow). Michael Villeret, et. al, "A New Digital Technique for Implementation of Any Continuous PCM Companding Law,", IEEE Int. Conf. on Communications, 1973, vol. 1, pp. 11.12-11.17. MIL-STD-188-113, "Interoperability and Performance Standards for Analog-to-Digital Conversion Techniques," 17 February 1987. "TI Digital Signal Processing Applications with the TMS320 Family", pp. 169-198. [From Joe Campbell; Craig Reese, cfreese@super.org; Sepehr Mehrabanzad, sepehr@falstaff.dev.cdx.mot.com] ============================================================================= Q2.7: How can I do CD <-> DAT sample rate conversion? CD players use a 44.1 kHz sample rate, whereas DAT uses a 48 kHz sample rate. This means that you must do sample rate conversion before you can get data from a CD player directly into a DAT deck. [From Ed Hall, edhall@rand.org:] For a start, look at "Multirate Digital Signal Processing" by Crochiere and Rabiner (see FAQ section 1.1). Almost any technique for producing good digital low-pass filters will be adaptable to sample-rate conversion. 44.1:48 and vice-versa is pretty hairy, though, because the lowest whole-number ratio is 147:160. To do all that in one go would require a FIR with thousands of coefficients, of which only 1/147th or 1/160th are used for each sample--the real problem is memory, not CPU for most DSP chips. You could chain several interpolators and decimators, as suggested by factoring the ratio into 3*7*7:2*2*2*2*2*5. This adds complexity, but reduces the number of coefficients required by a considerable amount. In any case, your local DSP guru will probably be able to suggest a good program for designing digital filters (this is a good general question for the net, anyway). This is one case where going beyond the simpler windowed-sinc-based approaches might pay off (and I'm afraid I can't recommend any readily available programs myself). ============================================================================= 3. Programmable DSP chips and software. This section deals with programmable DSP chips and their software. Q3.1: What are the available DSP chips and chip architectures? { This is based on a woefully inadequate databook collection. Anyone want to add to this list? Manufacturers want to submit anything? } The "big four" programmable DSP chip manufacturers are Texas Instruments, with the TMS320 series of chips; Motorola, with the DSP56000 and DSP96000 series; AT&T, with the DSP16 and DSP32 series; and Analog Devices, with the ADSP2100 series. A good overview of prorammable DSP chips is published periodically in EDN magazine. The most recent version is from Sep. 17, 1992, pp. 90-141. Here's a less ambitious chip breakdown by manufacturer: ------------------------------------- Texas Instuments: TMS320C1x: family of low cost fixed-point DSP's; 16 bit data, 32 bit registers; Various RAM and ROM configurations; 16 bit I/O bus, serial ports. TMS320C25: 50MHz fixed-point DSP; 16 bit data, 32 bit registers; 12.5 MIPS @ 50MHz. TMS320C30: 27/33/40 MHz floating point DSP; 32 bit floating point, 24 bit fixed point data, 40 bit registers; DMA controller; dual serial ports; some support for multi-processor arrays. TMS320C31: version of C30 minus peripheral bus, one serial port, and the 4Kx32 internal ROM. ~$20, 132 pin PQFP. TMS320C40: 40/50 MHz floating point DSP; extensive parallel processing support through 6 buffered byte-wide 20 Mb/s links and 6 channel DMA; cache. TMS320C50: enhanced TMS320C25 (double throughput); low overhead looping; 10 Kwords SRAM on chip. ------------------------------------- Motorola: DSP56001: 20.5, 27, or 32 MHz fixed point DSP. 24 bit data bus, 16 bit address bus, 56 bit accumulators (2), host interface port, serial ports (2), general purpose I/O pins, timer. Harvard architecture. 512 words program RAM, 512 words data RAM on chip. DSP56000: Mask-programmed version of DSP56001. DSP56002: DSP56001 with On-Chip Emulation (OnCE) debug port and clock PLL. Also has a four cycle double precision multiply and support for block floating point. Available up to 40 MHz. DSP56156: 40, 50, or 60 MHz fixed point DSP; 16 bit data bus, 40 bit accumulators (2), host interface port, serial ports (2), timer, OnCE debug port, clock PLL, 14 bit sigma-delta voice band CODEC, 2K words program RAM, 2K words data RAM on chip. DSP96002: IEEE format floating point DSP; 32 bit data and address bus, two complete external buses; Harvard architecture. ------------------------------------- AT&T: DSP32C: floating point DSP; 32 bit floating point, 16/24 bit fixed point data. DSP3210: floating point DSP; 32 bit floating point, 16/32 bit fixed point data; 32 bit address and data bus, serial port. [and others; we don't have a good list, unfortunately]. ------------------------------------- Analog Devices: ADSP2100: 32 and 50Mhz fixed point DSP (8 MIPS, 12.5MIPS). 16 bit registers except for multiplyer-accumulate register which is 40 bits. No on chip memory except for a 16 word instruction cache. ADSP2101: Derived from ADSP2100, 16 bit registers except for the multiplier accumulator which is 40 bits. 2Kx24 instruction/data ram in program memory space, 1Kx16 data ram in data memory space. Adds memory, timer, serial ports, etc. to the 2100. Fastest speed grade in production is 16.6 MHz (16.6MIPS). ADSP2102: Ram/rom version of 2101; user selects how much of the 2kx24 program memory is mask rom. ADSP2103: 3V version of the 2101. ADSP2105: 10Mhz fixed point DSP with 1 serial port, timer and 1kx24 instruction/data ram in program memory space, and 512 word data ram in data memory space. This processor sells for US $9.90 in any quantity. ADSP2111: adds a 8/16bit host interface port [to 2101?]. ADSP21msp50: ADSP2111 with an on chip a/d and d/a interface and additional low power modes. ADSP2161: has 8Kx24 mask rom. DSP21020: 20/25/33 MHz floating-point DSP; Supports 32-bit fixed point, IEEE format 32-bit floating point, and 40-bit floating point; 40-bit registers plus two 80-bit fixed-point multiply-accumulators; Harvard arch. with 32 word instruction cache allows two data accesses in a single cycle; IEEE 1149.1 JTAG boundry scan; 33.3 MIPS @ 33.3 MHz. ADSP21010: Slower and cheaper version of '020 (16 MHz). Limited to 32-bit fixed and floating point. All of the processors (except the 2100) use a 1X instruction clock and use an on chip PLL to generate an internal 4X clock. All processors have an extended Harvard architecture which allows two data fetches and an instruction fetch every cycle in parallel with an alu or mac operation. All instructions including accessing external memory can complete in 1 cycle. [Greg Koker, greg.koker@analog.com] =======================================================================